Highly controllable qubit-bath coupling based on a sequence of resonators 
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Combating the detrimental effects of noise remains a major cliallenge in realizing a scalable 
quantum computer. To help to address this challenge, we introduce a model realizing a controllable 
qubit-bath coupling using a sequence of LC resonators. The operating principle is similar to that 
of a recently proposed coplanar-waveguide cavity (CPW) system, for which our work introduces 
a complementary and convenient experimental realization. The lumped-element model utilized 
here provides an easily accessible theoretical description. We present analytical solutions for some 
experimentally feasible parameter regimes and study the control mechanism. Finally, we introduce 
a mapping between our model and the recent CPW system. 
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I. INTRODUCTION 

Achieving a balance between the seemingly compet- 
ing demands of qubit addressibility and coherence time 
remains a fundamental obstacle on the path to engineer- 
ing a large scale quantum computer[l|-[3|. By capitaliz- 
ing on the fundamental properties of photon-matter in- 
teraction in the solid state, the field of circuit quantum 
electrodynamics (cQED)U^ provides an attractive ap- 
proach to removing this obstacle. Recent advances in the 
field have led to the realization of quantum bits, qubits, 
approaching the requirements set by quantum threshold 
theorems [lO|, lll| and shown that coupling a qubit to a 
coplanar-waveguide (CPW) resonator can offer protec- 
tion against unwanted noise |12Hl4| . A promising demon- 
stration of the control over the system-bath interaction 
provided by the cQED architecture is the proposal for 
studying single-photon heat conduction(l5l. Il6l| . Entering 
a regime where such a conduction mechanism is dominant 
requires accurate tuning of the system, accompanied by 
extreme robustness against noise. 

A recent proposal introduced a cQED setup in which 
the coupling between a superconducting qubit and a ther- 
mal environment may be tuned at will, allowing for ei- 
ther rapid initialization or protected evolution of the 
gubitflTfl. Such a setup offers several advantages over 
previous works |18H20l | and is, in principle, realizable with 
current technologies. The setup is based on a CPW trans- 
mission line divided into two capacitively coupled cavi- 
ties, one housing a resistor and the other housing a qubit. 
By carefully selecting the system parameters, we can en- 
ter a regime where the resistor acts as the dominant noise 
source for the qubit and the coupling between the two 
can be externally tuned by manipulating the Josephson 
inductance of a set of superconducting quantum inter- 
ference devices (SQUIDs) placed inside the cavity. The 
proof of principle for this tunable environment exploited 



the distributed element (DE) model, wherein the CPW 
cavity is represented by an infinite number of differen- 
tial capacitors and inductors [21|. The model also allows 
one to account for the positioning of the resistor and the 
qubit inside the cavity. However, it would be beneficial 
to find a model system exhibiting similar tuning charac- 
teristics, but in which the complete dynamics are more 
easily accessible and some of the requirements of the DE 
model are lifted. 

In this paper, we study a setup where the CPW cavities 
are replaced by a pair of coupled quantum LC oscillators. 
In our analysis, and unlike in Ref. [13, the resistor is ca- 
pacitively rather than galvanically coupled to one of the 
oscillators which is an advantageous property in the fab- 
rication point of view of such devices. The result of this 
study is a simple and intuitive model for a qubit with an 
externally tunable coupling to a dissipative environment 
over many orders of magnitude. Furthermore, the model 
system has an analytical solution for several experimen- 
tally important cases. As a means of comparison, we es- 
tablish a meaningful mapping between the model system 
studied here and the DE model system, and demonstrate 
that we can arrive at the same eigenspectrum. 

The outline of the paper is as follows. In SecUll we in- 
troduce our physical system. Analytical solutions for the 
low-energy dynamics are given and analyzed in Sec. IIIII 
with special attention devoted to the protection of the 
qubit. In Sec. IIVI we introduce the operating principle 
of our system and study the control over the qubit life- 
time. Section |V] presents a mapping between the CPW 
cavity model system and our model system. We compare 
the excitation spectra between the models and study the 
emerging effective resonator-qubit coupling capacitance. 
We conclude the paper in Sec. I VII 



II. SYSTEM 
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Our system includes a sequence of lumped-element LC 
resonators. For simplicity, we restrict our analysis here to 
two such resonators and denote them as left and right cor- 




FIG. 1. Effective circuit diagram of our model system. Tlie 
left and right LC resonators are coupled by a capacitor Cc- 
A resistor of resistance R is capacitively coupled to the left 
resonator via a capacitance Ce and acts as a noise source. 
The right resonator is coupled via capacitance Cq to a Cooper 
pair box whose Josephson capacitance is Cj and Josephson 
energy Ej. 



tween them. The resonator Hamiltonians obtain an ef- 
fective form due to the internal system couplings 



Hk 



CkVl: + -Lull 



k = L,R, 



(2) 



where Vfe is the operator for the voltage across the kth 
LC circuit, with j^. the corresponding current operator. 
The effective capacitances are Cl = Cl + Cc + Ce and 



Cr — Cr 
term is 



Cc + Cq. The resonator-resonator coupling 



Hr 



L-R 



-CcVlVr, 



(3) 



and, similarly, the coupling between the resonators and 
the resistor assumes the form 



i?i, 



-CeVlSV,, 



(4) 



responding to the schematic representation of the circuit 
given in Fig. [T] Each resonator is comprised of a lumped 
inductor of inductance L^ and a lumped capacitor of ca- 
pacitance Cfc, where k = L,R refers to the left and right 
resonators. The resonators are coupled through a capac- 
itance Cc- To introduce a coupling to a bosonic environ- 
ment, a resistor of resistance R is capacitively coupled to 
the left resonator via the capacitance Ce- The right res- 
onator is coupled to a Cooper pair box (CPB) operated 
as a qubit in the transmon regime [12J. This coupling ca- 
pacitance is denoted by Cq, and the Josephson junction 
defining the superconducting island has capacitance Cj 
and Josephson energy Ej- In experimental realizations of 
the system, the inductance of the left resonator is made 
tunable by inserting one or more SQUIDs, whose Joseph- 
son inductances are tuned by external magnetic fields, in 
series with the inductor. 

On the general level, the working principle of our sys- 
tem described in Fig. [T] is the following. The left res- 
onator acts as a tunable filter for the noise and dissipa- 
tion arising from the resistor. Thus if the left resonator 
is tuned in resonance with the qubit, the noise at this 
frequency is transmitted through the left resonator and 
introduces dissipation for the qubit. If the left resonator 
is in resonance with the right resonator, strong dissipa- 
tion is introduced there. If the left resonator is far off 
resonance with either the right resonator or the qubit, 
the noise is heavily filtered and the effect of the resis- 
tor on the dynamics of the corresponding element is very 
weak. 

The Hamiltonian for our model can be written as the 
sum of the charging, inductive, resistive, and Josephson 
Hamiltonians corresponding to each element in Fig. [TJ If 
we neglect terms of second order in the resistor voltage 
fluctuations, the Hamiltonian becomes [22| 

Htot = HL+HR + HL^R+Hint+Hrcs+Hq + HR^q, (1) 

where each term describes a subsystem or a coupling be- 



where SVrcs describes the intrinsic voltage fluctuations 
across the resistor [22|. The resistor Hamiltonian, ffres, 
represents a bosonic thermal bath [23. The resistor can 
also be modelled as a transmission line [2J| although this 
approach is not employed here. The qubit Hamiltonian 
consists of the charging and Josephson parts (25| 



Hq = ^{Cj + Cq)V^^+Ej cose 



(5) 



where Vq and (/) are the operators for the voltage and 
superconducting phase difference across the junction, and 
for simplicity, we have neglected the self-capacitance of 
the superconducting island. The coupling between the 
right resonator and the qubit is given by 



Hr. 



-CgVRVq. 



(6) 



The resonator Hamiltonians given in Eq. ^ are of the 
form of quantum harmonic oscillators and are diagonal- 
ized by introducing a J, {cik) as the relevant bosonic cre- 
ation (annihilation) operator for the kth resonator. The 
diagonalization yields 



t;^ 



Hk = hujk{a\ak + 1/2) 



(7) 



where ujk = l/-\/LfcCfe, and allows us to define the voltage 

operators as Vk = V^?(a[.-l-afe), where VjP = \JfujJk/{2Ck): 

and the current operators as Ik = /^(aj, — Ofc), where I^ = 

iy^hajk/{2Lk)- Using the same formalism, the resonator- 
resonator coupling term can be written as 



^t;;t 



Hl-r = ha{a\a\ + a\aR + Sllci^ + cllcir), (8) 

where we define a = —CcVj^V^/ha.s the resonator- 
resonator coupling strength. Assuming that the 
qubit is operated in the transmon regime, the op- 
erator for the voltage across the junction is Vq ~ 

TTk i^) ^" (^ - ^^)' ^ith Ec = eV[2(C^ + Cq)\ 



'Cj+Cg \SEc J 



being the charging energy for the CPB island and ¥ and 
b being the creation and annihilation operators for the 
harmonic oscillator approximating the transmon, respec- 
tively |12| . Hence, the resonator-qubit coupling Hamil- 
tonian becomes 



excited states are 



H 



R-q 



hg{ba\ + ban - V^an - ^'''al), 



(9) 



. V2£C^ / E,, \ 
''h{Cj+Cg) \SEc J 



1/4 



is the resonator-qubit 



where g 
coupling strength. 

Our main interest is the dynamics of the resonator- 
qubit system and, to this end, we treat the resistor as 
a dissipative environment described by the thermal bath 
Hamiltonian iJres- Assuming that the coupling to the 
environment is weak, the coupling term iJjnt can be taken 
as a small perturbation inducing transitions between the 
resonator-qubit states defined by the eigenproblem (H^ + 
Hjl + Hl^h + Hq + Hii^q)\rn) = E,ri\rri). The transition 
rates are approximated using the Fermi golden rule|26| 
as 



|(n|CgVL|m)p 



■SsV,^S~^rnn)- 



(10) 



Here, ujmn = (£'« — Em)/h and we assume that the resis- 
tor is in thermal equilibrium so that the voltage fluctu- 
ations (5VI-CS follow the Johnson-Nyquist spectral density 
Ssv.^^oj) = 2hRuj/[l - e-'i'^/(fcB^)], where T is the effec- 
tive resistor temperature. 



III. ANALYTICAL RESULTS 

Attaining closed-form solutions for the eigenproblem 
presented in the previous section is likely not plausible 
in the general case and we therefore focus our efforts on 
a few important special cases in the low-energy regime. 
This analysis provides convenient tools for experimen- 
talists working on the system and offers insight on the 
robustness of the qubit against noise. 



A. Decoupled qubit 

We begin by studying the dynamics of the system 
shown in Fig. [T] when the qubit has been completely de- 
coupled, i.e., when Cg = 0. We define the excitation 
number product state basis as {|ni, rifl)}, where n^ and 
ur are the excitation numbers of the left and right res- 
onators, respectively. Furthermore, we consider only the 
low-energy dynamics so that the basis is restricted to 
{|0,0), |0, 1), |1,0)}, omitting simultaneous excitations of 
the left and right resonators. We choose the zero-point 
energy so that the lowest eigenenergy is i?o = corre- 
sponding to the ground state \g) — |0,0). The subse- 
quent excited state energies are E± ~ fi,(f7± VA^ -f a^), 
where fi = (wi -I- ujji)/2 and A — [u^ — ujr)/2^ and the 



I-) = -cos(6l)|l,0)+sin(6l)|0,l), 
1+) =sin(6')|l,0)+cos(6l)|0,l), 



(11) 
(12) 



where tan(6') = (A + VA^ +a^)/a. With the help of 
Eqs. pT?)l -p ^ . the relaxation rates can be written as 



r__ 



a 



Ruj- 



Ce_ 

[(A + VA2 + a2)2 + Q,2] Zl [1 - e-'^'^-/(kBT^ y Cl 

(13) 

(A + VA2+a2)2 J?a;+ (Ce 

^^' [(A + VA2 + a2)2 + ^2j Zl [1 - e-''-+/(^-^)] ^ Cl 

(14) 

and the corresponding excitation rates are acquired 
through a substitution Vg^±{uj) = V±^g{—uj). Above, 
we defined the characteristic impedance of the fcth res- 
onator as Zk = yLkJCk and used a notation lu± ~ 
E±/h. Near resonance in which A — > 0, we have sin(6') ^• 

l/-s/2 and cos(0) — > l/v^, and hence r+_s.g « F yg 

for a <C ri. The rates at the resonance are Flfij. = 

R{n ± a){CE/CL)y[2ZL{l ~ e-'^("±")/('=«^))] and the 
eigenstates become equal superpositions of the excited 
resonator states. 

In the far-detuned limit A 2> a, we have cos2(6') « 



a2/(4A2) and sin^(6') w 1 - a^/{4A^), in the case of 



which the transition rates reduce to 



a 



Rlo^ 



-^9 ~ 4^2 Zl[1 - e-fi--/(feBT)] ycL 

^2 X ^^^ 



Ce 



nfar 



1- 



a 

4A^ 



ZLil-e-^^+ZC^BT)^ \C_ 



Ce 



(15) 



(16) 



The eigenstate |— ) approximately corresponds to an ex- 
citation in the right resonator whereas the eigenstate | +) 
corresponds to an excitation in the left resonator. To 
quantify the effect of the detuning, we can define the ratio 
between the decay rate of an excited state in Eqs. ((H])- 
(|12p far from resonance and at resonance, the tuning 
ratio r± = T^^'i^g/T'i^'L^g. For the state |— ) approxi- 
mately corresponding to the excited state of the right 
resonator far from resonance, the ratio is approximately 
r_ « a2/(2A2) for a, A <C f^. Hence we can effectively 
decouple the right resonator from the noise source by de- 
tuning. As a consequence, a qubit directly coupled only 
to the right resonator can be efficiently protected. 



B. Resonator-qubit system 

We proceed by analyzing a few important regimes 
where a closed-form solution for the full system in Fig. [1] 
is accessible. This analysis not only sets the premise for 
studying the general behavior of the system, but also of- 



fers simple tools for experimentalists operating in these 
important regimes. We extend the computational ba- 
sis from the previous section to include the qubit so 
that the basis states are of the form \nL,nji^a), where 
a = g,e for the ground and excited qubit states, respec- 
tively. The basis relevant to the low-energy dynamics be- 
comes {|0, 0,5), |1,0, (7), |0, l,.g), |0,0, e)} which does not 
include states of the form of tensor products of the ex- 
cited qubit state and resonator excitations as those only 
couple to higher-order photonic excitations. This omis- 
sion invokes effectively the rotating wave approximation 
so that the resonator-qubit interaction is of the form 
Hji^q = h{ga]^a~ + g*aiia^), where (t+ and a~ are the 
Pauli raising and lowering operators for the qubit, respec- 
tively. Within the selected basis, the time-independent 
Schrodinger equation for the resonator-qubit system can 
be written in a practical form yielding both the eigenen- 
ergies and eigenvectors [see Appendix R]. 

First, we study the partial resonance case where ujr 7^ 
lul ^ iOq = uj, that is, the qubit is in resonance with the 
dissipative resonator but not with the one to which it is 
directly coupled. Here hiOq denotes the qubit excitation 
energy obtained by diagonalizing Hq. Using Appendix El 
the excited state energies {efe} obtain closed- form solu- 
tions as ei = hjj, £2/3 = fi-(wav± \/^'d + l.9p + "^^)' where 
Wav = ('^ + '^_r)/2 and Wd = (w — W7j)/2. The correspond- 
ing eigenvectors are 



\Vi) 



\'f2) = 



2 \a 



"+ Iff 



(||l,0,5) + |0,0,e)), (17) 



A+l"^ 



m) 



M|^,,„,.„-. 






|0,l,.g) + |0,0,e) 
(18) 

|0,l,.g) + |0,0,e) 
(19) 



where A = ^j {ud/\9\f + {a/\g\f + 1 and Vt± = 

\J a^ + IgP -I- i^A\g\ ± ujd)^ . Note that the ground state 
is lyjo) = lOiOiff) with the zero-point in energy chosen so 
that eo = 0. Using Eqs. P?)l - (|19p . the relaxation rates to 
the ground state become 



T^qrcs 



RiO\ 



|5|2 -f a2 Zi\\ - e-'i'^i/(feBT)] \^Cl 



"¥] . (20) 



T^qrcs 



a 



RixJo 



Cf 



nl Zl[1 - e-f"^^-/^''^'^')] \c_ 



pqr 



a 



Rus, 



^^^"^f n\ Zl [I - e-'*-'3/(fcBT)] \^ Cl 



Ce 



(21) 
(22) 



where 



^ml'h-, m = 1,2,3, for simplicity. 



The excitation rates are obtained from F^^Jf^.^^ (w) = 
ry'L,^ (— w)- The full resonance solution is obtained 
by setting oj^ ~ uj and if we additionally set a ~ 



|g|, the rates satisfy V^^^^ 



and obtain the form F^®^% 



when v2a ^ u, 

Ru:{CeICl?/[^Zl{1 - e-'^-/^'^--^))] and F^^^^„ _ 
i?(a;±v^a)(C£;/CL)V[4ZL(l-e-'*("±^")/('=«^))] with 
plus and minus corresponding to /c = 2 and fc = 3, re- 
spectively. 

If we study the case where the resonator-resonator cou- 
pling is much stronger than the resonator-qubit coupling, 
i.e., a ^ \g\. the closed-form solutions indicate that the 
qubit decay is much slower than the decay of the res- 
onators. Here we identified |(/3i) as the excited state of the 
qubit and \ip2) and |(/?3) as superpositions of the excited 
states of the resonators by studying the corresponding 
state amplitudes. If a ^ jgl, on the other hand, the de- 
cay of the left resonator can be seen dominate over both 
the qubit decay and the decay of the right resonator. 
Similarly to the previous case, we identified \ipi) as the 
excited state of the left resonator and 1932) and \ip^) as 
superpositions of the excited states of the right resonator 
and the qubit. Note that these observations in the res- 
onance case serve as a consistency check for the system 
and are not to be confused with the control introduced 
by manipulating the resonator-resonator and resonator- 
qubit detunings. If a and \g\ are of the same order of 
magnitude, the detunings determine the relative decay 
rates as we will show in the following section. 

Next, we study the case where the system couplings 
are weak \a\, \g\ ^ cj^, cjj^, ujq and the system is far away 
from resonance |a|, \g\ <C jw; — Ws|, for all Z, s S {L, R, g}, 
so that a perturbation theory solution becomes available. 
After expanding to second order in the small coupling pa- 
rameters, the resulting eigenenergies and states are given 
' in Appendix [Bl Using these results, the relaxation to the 
ground state takes place at rates defined by 



F 



far 



- 42 

- ^1 



/2A 



LR 



\ 2A2 



LR 



Ru, 



Ce 

(23) 



-ifar 



^Ai 



Ru!o 



Ce 



pfar ^ 42 ISP"^ 



i?(x)3 



^lR^lLZL[l-e-'-^^^n^BT)^^ \Cl 



Ce 



(24) 



(25) 



Vs-^Vo^ 



Here we denote Ajs = w; — w^ as the difference between 
two intrinsic system frequencies and the normalization 
constants A^ are given in Eqs. (jB7p - (|B9[) . Note that 
Ak = l + 0(|a|/|Ais|,|5|/|A/,|), for all k =1,2,3 and 
l,s = L,R,q. Assuming |a|,|5| <C |w; — w^^j implies 
that Ak ^ 1 making Eqs. ((23)) - (p5| easier to interpret. 
We observe that the relaxation rate of the state {(p^) 
(« |0, 0,e)) is proportional to [gpa^, that is, the excited 
qubit state is doubly protected against noise: firstly by 
the weak resonator-qubit coupling and secondly by the 
weak resonator-resonator coupling. This is in contrast to 
state \(p2) (~ 10,1,(7)) which is proportional to a^ and, 
hence, only benefits from the weak resonator-resonator 



coupling. Similarly to Sec. IIII Al we can define the tun- 
ing ratio r^3 = T^^^L,^^ /r;^;%^^ for the state {ips) ap- 
proximating the excited state of the qubit far away from 
partial resonance. Using Eq. ^^ for the rate at partial 
resonance and Eq. (|25p for the rate far away from res- 
onance, we obtain r^g 
weak system couplings. 



|.g|VA^i for |wd| < |wav| and 



IV. TUNING THE QUBIT AND PHOTON 
LIFETIMES 

The control over the qubit lifetime stems from the 
resonators intermediating the relevant interaction with 
the bath. This becomes clear when we assume that the 
inductance of the left resonator can be externally con- 
trolled in the physical realization of the system. For in- 
stance, if the resonator describes a CPW cavity as dis- 
cussed in Sec. |Vl such manipulation can be implemented 
by the insertion of SQUIDs into the cavity. Manipulat- 
ing the frequency of the left resonator in this manner 
effectively allows us to drag the left resonator on and 
off resonance, lol ~ ^q being the resonance condition, 
with the qubit. A simple qualitative explanation for the 
resulting changes in the qubit lifetime can be found by 
considering the impedance of the LC resonators in a clas- 
sical circuit [231 . Far above resonance, the impedances of 
the resonators to the ground may be considered to be 
purely capacitive. Consequently the coupling-capacitor- 
resonator circuit acts as a voltage divider, resulting in 
effective voltage fluctuations which are reduced at the 
position of the qubit protecting it from decay. On reso- 
nance, the impedance of the resonator to the ground may 
be considered infinite so that the qubit-bath coupling is 
of the form of an effective series capacitance resulting in 
fast decay. 

In the following simulation, we identify two of the low- 
energy eigenstates of the qubit-resonator system as the 
ones approximating \q) « |0, 0, e) far from the resonator- 
qubit resonance ujq = wl, and \R) « |0, 1, g) far from the 
resonator-resonator resonance un = u^. Even though 
these eigenstates do not match the excited states of the 
qubit and the right resonator near the resonances, they 
still allow us to study the operating principle of the sys- 
tem. We calculate the decay rate of \q) to the ground 
state, Fg, from Eq. pOp by numerically diagonalizing the 
Hamiltonian in Eq. ([T]), and present the rate in Fig.[5]as a 
function of the inductance of the left resonator. The rate 
obtains a sharp peak when the left resonator is tuned 
near resonance with the qubit, as expected. Note that 
with the realistic system parameters used in Fig. [51 tun- 
ing the inductance allows for significant control over the 
qubit lifetime. For instance, adjusting the resonance fre- 
quency of the left resonator from the point lol = oj^'^'^^, at 
which Tq obtains its maximum value, to uj^/oj^'^^'^ = 0.75 
yields a decrease in the decay rate given by the ratio 
F,(wP<=^'^)/Fg(0.75wP°^'') = 1.69 X 10^. Hence, externally 
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FIG. 2. (Color online) Relaxation rate to the ground state 
of the eigenstate |g), F,, and the eigenstate |_R), Vr, as a 
function of the inductance of the left resonator presented by 
the solid lines. Far from resonance, these rates are accurately 
approximated by the weak-coupling solutions of Eqs. H24p - 
(|25p presented by the dotted lines. For the coupling capaci- 



capacitors of magnitude Cl ~ Cr = 39 pF, and the induc- 
tance of the right resonator is fixed to Lr — 40 nH resulting 
in u}r/{2-k) = 12.79 GHz. The resistor has a resistance of 
R — 500 Q, and an effective temperature of T = 10 mK, and 
the qubit parameters are chosen such that Ej/Ec = 50, while 
its frequency is fixed at Uq/{2Ti) = 15.92 GHz. The resonator- 
qubit coupling strength is |(?|/(27r) = 178.25 MHz and the 
resonator- resonator coupling strength ranges from Q/(27r) = 
-19.40 MHz at Ll/Lr = 0.5 to a/(27r) = -14.74 MHz at 
Ll/Lr = 1.5. 



rapidly initialize the qubit to the ground state and to 
protect it from dissipation. The frequency ojP'^^'*^ differs 
slightly from the qubit frequency ojq due to the resonator- 
resonator coupling strength having a nonlinear depen- 
dence on lol. The analytical solutions for when the sys- 
tem couplings are weak given in Eqs. ([23)) - p5|) accurately 
approximate the numerical solutions far away from reso- 
nance in Fig. [21 In particular, the rate F^ is accurately 
approximated by F^^^^ in Eq. ([55)1 near the resonator- 
resonator resonance despite the formal assumption that 
the system remains far away from all resonances used 
when deriving the analytical solutions. Even though non- 
idealities in the experimental realization of the system as 
well as other dissipation mechanisms can limit the phys- 
ical rates in practice, the decay caused by the coupling 
to the thermal bath can be efficiently controlled in this 
manner. 

We show also the decay rate of \R) to the ground 
state, Vr. in Fig. [51 Similarly to the qubit case, this 
decay rate exhibits a strong peak which, however, re- 



sides at lol 



peak 



'R 



near the resonator-resonator res- 
onance ujl = ujR. Comparing the decay rates us- 
ing the same relative changes in the frequency of the 

1 -^ 0.75) yields a ratio of 



left resonator {ojl/oj^'^ 
F«(c.r')/rfl(0.75wr') 



2.66 X ID* where 



^L 



peak 



manipulating the inductance makes it possible to both is the frequency at which Tn reaches its maximum value. 



Even though the range of control is more hmited than 
in the qubit case, we can still adjust the decay rate 
over several orders of magnitude. The difference in the 
ranges is an indication of the double protection of the 
qubit against noise mentioned in Sec. lIIIBl According to 
Eqs. ((24 |) -([25 | . the change in both Tq and Tr is approx- 
imately proportional to a^ as the inductance of the left 
resonator is manipulated far from resonance. However, 
Tq has an additional constant prefactor |gp allowing for 
stronger protection and, hence, a larger range. By set- 
ting Ljq = lul = w/j, we can rapidly initialise the entire 
system. Alternatively, by detuning both the resonators, 
as well as the qubit, we can protect both the excited right 
resonator state, and the qubit from decay. 



MAPPING FROM THE DISTRIBUTED 
ELEMENT MODEL 



In Ref. [13, a similar-purpose physical system is consid- 
ered as here but composing of a coplanar- waveguide cav- 
ity and a qubit. The cavity is separated into two distinct 
regions by inserting a capacitor of capacitance Cc into the 
central conducting strip. These regions then house the 
resistor on the left and the qubit on the right along with 
an array of SQUIDs used to manipulate the inductance 
of the left region. The application of the distributed ele- 
ment (DE) model to this situation allows one to account 
for the positioning of the resistor and the qubit anywhere 
in the divided cavity. The key difference between this 
system and the one presented here is that the resistor is 
galvanically, as opposed to capacitively, coupled to the 
transmission line and therefore couples through the cur- 
rent rather than the voltage of the cavity. However, the 
setup using the DE model can be modified to incorpo- 
rate a capacitive bath coupling physically justifying the 
search for a meaningful mapping between the models. 
We present a simplified schematic of the cavity system 
with a capacitive bath coupling in Fig. |31 

The task at hand is to select the parameters such that 
the two systems accurately exhibit similar physical prop- 
erties. In the following, we refer to the resonators describ- 
ing the individual cavity regions as bare cavities. In order 
for the DE and LE models to yield the same spectrum for 
the isolated resonator and cavity systems, we match the 
frequencies of the lowest bare cavity modes to the res- 
onator frequencies in our model, and the corresponding 
voltage operators for the bare modes at the ends of the 
cavity regions to the resonator voltage operators in our 
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FIG. 3. Schematic representation of the CPW cavity system. 
The center conductor of the cavity divided by the coupling 
capacitor Cc placed at x — Xc is represented by the shaded 
elements. For clarity, the ground plane surrounding the cen- 
ter conductor is not shown. The capacitive bath coupling is 
introduced by coupling the left region via Ce to the resis- 
tor R. The right cavity region is coupled via Cq to a CPB 
defining a qubit with bare excitation frequency cuq. The qubit 
location x = Xq represents the position at which the qubit is 
connected between the center conductor and the ground plane 
of the cavity. 



model [see Appendix [C]. This yields a mapping given by 
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(26) 
(27) 
(28) 
(29) 



where c is the constant cavity capacitance per unit 
length, II {(-b) is the left (right) cavity inductance per 
unit length, Xc is the position of the dividing capaci- 
tor Cc in the central conducting strip of the cavity and 
Xcav is the length of the cavity. Including either the 
dividing capacitor in the CPW cavity or a resonator- 
resonator coupling capacitor of equal capacitance into 
our system produces matching low-energy spectra. An 
example of tuning the spectrum by manipulating the left 
cavity inductance per unit length is given in Fig. U) For 
the cavity parameters used in Fig. 21 the bare cavity 
resonance uj\ = o;^ is achieved at ^l = ^r [see Ap- 
pendix [C]. Away from resonance, the spectrum is that 
of isolated bare cavities. Near resonance, the capacitive 
coupling of the cavity regions induces an avoided crossing 
in the spectrum near which the bare cavity eigenstates 
are mixed. The minimum energy gap Ae™™ has a value 
of Ae'"'"/(27rft) = 32.82 MHz. 

Due to the fact that we match the voltage operators 
at the ends of the cavity regions to our resonator voltage 



operators, placing the qubit at Xq 



inside the right 



cavity region will give matching spectra for the cavity- 
qubit system and our resonator-qubit system. Note that 
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FIG. 4. (Color online) Excited state eigenfrequencies Wfc = 
ek/h, k = 1,2 (from bottom to top), as a function of the 
left cavity inductance per unit length. The eigenenergies e^ 
are calculated for the cavity system ignoring both the resis- 
tor and the qubit, and they are equal to those calculated for 
our system using the mapping detailed in this section. The 
eigenfrequencies are scaled by the frequency of the bare right 
cavity to'^ = 7r/[(a;cav — Xc)V£rc\. We use Xcav = 12 mm, 
Xc = Xcav/S and Cc = 1 fF. The cavity has a capacitance per 
unit length of c = 130 pF/m, and the right cavity region has 
a characteristic impedance of Zc = \JTeJc = 50 f2 and a bare 
frequency of J'rI{2-k) = 12.82 GHz. 



FIG. 5. (Color online) Effective resonator-qubit coupling ca- 
pacitance in the LE model with respect to the position of the 
qubit in the right cavity region in the DE model chosen such 
that the low-energy excitation spectra of the two models are 
equal. The solid line shows the numerical solution obtained 
through matching the low-energy spectra of the LE and DE 
models. The dashed line gives the approximative effective ca- 
pacitance C\ = Cq cos[7r(a::q — Xc3.v)/{xca.v — Xc)]- We use the 
same parameters as in Fig.|4]with a fixed II = £r- The cavity- 
qubit coupling capacitance in the DE model is Cq = 10 fF, 
and the qubit excitation frequency is ujq/{2n) — 12.62 GHz. 



the DE model includes the qubit only after the divided 
cavity eigenmodes have been solved and, hence, the ef- 
fective left capacitance Cl in the LE model must not 
include Cq in order for the two models to yield the same 
dynamics [2^. If an offset in the qubit position is de- 
sirable, we may define an effective capacitance C° in 
the LE model, selected such that the spectra of the DE 
and LE models still match. That is, we choose C^^ such 

that the eigenproblems {Hca,v + Hc-q + Hq)\Lp) = e^p\Lp), 
where Hca.v is the DE cavity Hamiltonian and Hc-q is the 
DE cavity-qubit coupling Hamiltonian, and [Hl + Hr 4- 
HL_ji+Hq+Hii_q)\ip') = e,p'|<p') give solutions for which 
{^f} = {^v'}- Figure [5] presents the dependence of C°* 
on the position of the qubit in the right cavity region. 
For Xq ~ Xcav, the effective coupling capacitance may be 



approximated by C°* = C° cos[7r(a;g 



.)/{Xc 



= )], 



where Cj? is the DE model cavity-qubit coupling capac- 
itance [29|. As shown in Fig. [51 the approximation is 
accurate in the vicinity of the exact matching point en- 
abling one to exploit the simpler LE model even if the 
qubit is not placed at the very end of the right cavity 
region. 



VI. CONCLUSIONS 

We have presented and analyzed a scheme for a su- 
perconducting qubit coupled to a bosonic thermal bath 
through a sequence of harmonic oscillator modes. The 
specific device studied in detail is based on a pair of cou- 



pled quantum LC oscillators intermediating the interac- 
tion between the qubit and the bath. We have shown 
that the qubit-bath coupling can be controlled by ma- 
nipulating the inductance of one of the resonators allow- 
ing for both rapid initialization to the ground state and 
protection from noise. 

Our design allows the use of lumped-element analy- 
sis rendering the anal ytic al results more appealing than 
previous descriptions |17| . It grants us access to simple 
analytic solutions in several experimentally interesting 
regimes. Furthermore, the setup we propose here makes 
use of a capacitive coupling to the thermal noise source, a 
feature which can be easier to fabricate than the systems 
in the previous proposals |15l - ll7| based on galvanic con- 
tacts. We propose a mapping between our scheme and a 
cavity-qubit design [l3| and find that the results from the 
two scenarios are in good agreement. Our model works 
as the proof of principle for obtaining tunability in qubit- 
bath interactions by implementing a coupling through a 
sequence of resonators. Even though our analysis concen- 
trated on a pair of LC resonators, we expect that a more 
general scheme in which the resistor is coupled to the 
qubit through N resonators offers increased tunability. 
Especially, if all the N resonators are out of resonance 
with each other and the qubit, we are likely to achieve 
iV-fold protection as opposed to the two-fold protection 
observed here in Eq. (PS)) . Finally, the analysis presented 
in this paper is not only directly usable for the experi- 
mental realization of such tunability but also for devising 
more elaborate designs based on similar ideas and appli- 
cations of these in various different physical systems. 
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Appendix A: Analytical solution to the energy 
eigenproblem 

In the low-energy regime and within the rotating-wave 
approximation, the time-independent Schrodinger equa- 
tion {Hl + Hr + Hl_r + Hg + HR_g)\<fk) = efc|<Pfc), 
where k £ {0,1,2,3}, has a ground state given by 
\ipa) = |0,0,g). See Eqs. (H])-® for definitions of the 
Hamiltonians. We select the zero-point energy such that 
eo — 0. After rewriting the above Hamiltonian in the 
basis described in Sec. IIIIBl the eigenproblem reduces to 
solving three simultaneous equations 



tan(7fc) 



tan(^fc) = -^ 
sin(7fc) 
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V\9?'{ei-io,y 



(Al) 
(A2) 

(A3) 



where ej. = Cfc/ft. The equations yield the eigenener- 
gies efc, k = 1,2,3, and, similarly, solving 7^ and 9k 
produces the excited states as \ipk) = sm{9k)\l,0,g) + 
cos{9k) sin(7fc)|0, l,.g) + cos{9k) cos(7fe)|0, 0, e). 



Appendix B: Weak-coupling limit 

Let us rewrite the total Hamiltonian of the resonator- 
qubit system [see Eqs. ©-(J!]) for definitions] ignoring 
the resistor as ffgys = ^l + ^_r + Hq + H', where 
H' = H^^f/ + Hj^^q is the small perturbation in the 
weak coupling limit |a|,|(7| ^ ujL,ujR,ujq and the sys- 
tem is far away from resonance |a|, Ig] ^ |w/ — w^l, for 
all l,s € {L,R,q\. The ground state is unaffected by 
the perturbation and given by \ipo) — |0, 0,(7) with its 
eigenenergy selected as the zero-point of energy. Af- 
ter applying the second-order perturbation theory in the 
small coupling parameters, the eigenenergies of i/gys are 
given by 



ei = hujL - 

e3 = hLJq- 



ha 






(Bl) 
(B2) 
(B3) 



where we adopted the notation A/^ = lui — ujs for the 
difference between two intrinsic system frequencies. The 
corresponding eigenstates are 



Ifi. 
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(B4) 

(B5) 

(B6) 



and the normalization constants are given by 



Ai^ 



A.= 



A.= 






ffl^Ai^ + 2a2|5PAi^A?,, + A^^^ {a^ + 4Ai^) 



2A^^A,i 



4A|JgPa2+4A4^^Ai^ + AiJg|^ 



(B7) 
(B8) 
(B9) 



Appendix C: Bare mode mapping 

If we omit the coupling capacitor inserted into the 
CPW cavity shown in Fig. [3l the cavity regions can be 
treated as individual resonators. We denote these res- 
onators as hare cavities and write the excitation frequen- 
cies for their lowest-energy modes as |2]| 



where (o;)' and o.;, I = L,R, are the creation and an- 
nihilation operators for the modes, respectively. The 

prefactors for the operators are y^° = \ fuJl/{xcc) and 



v^° = v^y[(2;c 



Xc)c]. 



{Xc 



(CI) 
(C2) 



where c is the constant cavity capacitance per unit 
length, l]^ {£fi) is the left (right) cavity inductance per 
unit length, Xc is the position of the dividing capacitor Cc 
in the central conducting strip of the cavity and a;cav is 
the length of the cavity. Similarly, the voltage operators 
for these modes are given by 



Vl{x)^Vl'co.^^y{ai)^ + a1] 



F^(x) = \/^"cos 



tt{x - 



[(4)^ + 4], 



(C3) 
(C4) 



We relate the bare cavities to our resonator system in 



Fig. [T] by demanding that uj^ 



ui 



Li 



^R = ^Ri ^1 



bO 
L 



yl 



and V^ = y^. This causes the resonators in our system 
and the bare cavities to have equal excitation energies, 
and the prefactors of the voltage operators at the ends 
of the bare cavities to match the prefactors of the volt- 
age operators in our system. The resulting mapping is 
given in Eqs. ([^5|) - (P9)l . Matching the voltage operators 
in this manner ensures that including the same coupling 
capacitors into both systems has the same effect on the 
spectrum. This is because the bare mode voltages at 
the coupling capacitor effectively determine its charging 
Hamiltonian. 
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